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It is argued that in heavy ion collisions at LHC there could be a sizeable production of baryons 
containing two or three heavy quarks from statistical coalescence. This production mechanism is 
peculiar of Quark Gluon Plasma and the predicted rates, in heavy ion collisions at LHC energy, 
exceed those from a purely hadronic scenario, particularly for H(,c and Hccc Thus, besides the 
interest in the discovery of these new states, enhanced ratios of these baryons over singly heavy 
flavoured hadrons, like B or D, in heavy ion collisions with respect to pp at the same energy, would 
be a clear indication of kinetic equilibration of heavy quarks in the Quark Gluon Plasma. 


Charm and bottom quarks are expected to be abun¬ 
dantly produced in hadronic collisions at very high en¬ 
ergies. In heavy ion collisions (HIC), multiple pair pro¬ 
duction is expected to occur, with average multiplicities 
which, at the LHC energy of 5.5 TeV, attain 0(10) for 
bottom and 0(100) for charm in central collisions. These 
quarks, produced in the early stage of the collision off 
hard scatterings, lose energy in the Quark Gluon Plasma 
(QGP) and, if the lifetime of the source is long enough, 
may reach thermal P| (not chemical, as their annihilation 
rate is very low) equilibrium within the medium. At the 
hadronization point, they will coalesce into hadrons. If 
coalescence process occurs statistically at the hadroniza¬ 
tion temperature, there is a finite chance that two or even 
three of them coalesce into the same particle, thus giv¬ 
ing rise to multiply heavy flavoured hadrons, particularly 
baryons. This phenomenon is likely to occur only if heavy 
quarks get very close to thermal equilibrium reshuffling 
over a large region because high momentum quarks, most 
likely, will hadronize into different particles unless two or 
three of them emerge very close in momentum from the 
hard process. The latter production mechanism, where 
multiply heavy flavoured hadrons arise from correlated 
quarks, predicts multiplicities which, at some large en¬ 
ergy, are exceeded by those predicted by coalescence of 
uncorrelated quarks. The ultimate reason of this effect 
is that the average multiplicity of heavy quark pairs in¬ 
creases faster than soft hadrons multiplicity as a function 
of centre-of-mass energy. This is in turn related to the 
volume of the system at freeze-out. Consequently, the 
system gets denser in heavy quarks at chemical freeze-out 
as energy increases and so does the chance of formation 
of multiply heavy flavoured hadrons. Therefore, an en¬ 
hanced production rate of these objects relative to singly 
heavy flavoured hadrons (like B’s or D’s) is distinctive 
of HIC and can be used as a probe of thermalization of 
heavy quarks within the QGP, thence as a signal of QGP 
itself. This idea was advocated in refs. [2, 0 where the 
authors envisaged an enhancement of Be and J/if} mesons 
production from the QGP with respect to hadronic colli¬ 
sions at RHIC. In this paper, we amend and reinforce this 
picture by observing that the main advantage of baryons 
(and especially and Qcce whose signal should be de¬ 


tectable at LHC) over quarkonia and Be is to enhance 
the difference from the “background” of coherent produc¬ 
tion, i.e. direct production of these states from early hard 
scatterings, possibly followed by melting in the plasma. 
Furthermore, in view of this relatively large production, 
HIC become a suitable place to discover many of these 
yet unobserved states. 

That the production of multiply charmed baryons 
could be enhanced in high energy HIC was first proposed 
in ref. 0 on the basis of a quark recombination model. 
A different coalescence picture was used in ref. . How¬ 
ever, in both studies, no relation was set between produc¬ 
tion enhancement and heavy quark kinetic equilibration 
in the plasma and quantitative predictions of production 
rates were dependent on an undetermined free parameter. 
On the contrary, predictions are definite within the statis¬ 
tical coalescence model (SCM), whose basic idea has been 
introduced in ref. 0 following a work on statistical pro¬ 
duction oi J/ij) and used thereafter by many authors 
0. In practice, the SCM is the statistical hadroniza¬ 
tion model supplemented with the constraint of a fixed 
number of heavy quarks and antiquarks. 

We are now going to calculate of the mean multi¬ 
plicities {uj) of heavy flavoured particles in the frame¬ 
work of the SCM. In HIC at high energy, electric charge, 
strangeness and baryon number conservation can be 
treated grand-canonically whereas charm and beauty 
conservation cannot because the multiplicities of heavy 
flavoured hadrons are not large. Thus, besides the num¬ 
ber of c -|- c quarks and A"t, of b -I- b quarks, also 

the net charm C and net beauty B should be fixed. In 
fact, instead of C, B, Nc, Nb, any combination of these 
integers can be used to constrain the partition function. 
Indeed, it is advantegeous to use the numbers Vc of c, Vc 
of c, Vb of b and vj, of b quarks. The relevant partition 
function reads: 


Z{iZc,Uc,Vb,Vb) = Zi 


n 


27r ' 




X exp 


f=c,c,b,b 
0c -ii'oj ipc-ii'bjipb-n'bj'Pb 


( 1 ) 






2 


where Zi is the grand-canonical partition function includ¬ 
ing all light-flavoured species; are the one-particle par¬ 
tition functions: 
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Xj are the fugacities (with regard to electric, baryonic and 
strangeness charge) and Vcji r'gj, vi^ are the number 
of c, c, b, b quarks respectively of the of the j**' hadronic 
species; the factor gj in Eq. Q is its spin degeneracy 
and T is the temperature. The sum over j in Eq. o in¬ 
volves all hadrons containing heavy quarks. The primary 
average multiplicity of any of these, in events with fixed 
numbers of heavy quarks, reads: 
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where: 
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and similarly for C{^f) with replacing the afn’s. We 
note in passing that = 0 if < 0 and C( 0 ) = 1 - 

The mass hyerarchy in heavy flavoured hadrons make 
a /1 a /2 a /3 (/ = c,b), because, according to lO, 

each term is suppressed with respect to the preceding 
one by a factor ~ exp(—m//r) where m/ is the c or b 
quark mass. Neglecting a /3 (i.e. flccc or flbbb baryons), 
the integrals in 0 give rise to a polynomial expression: 
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It can be realized that we would deal with factorized 
charm and bottom integrals in Eq. o were not for the 
presence of hadrons carrying both flavours, such as Be. 
To recover factorization, we can expand the integrand in 
Eq. in power series of the Zj of hadrons containing 
both c and b (anti-)quarks up to the first order. Defining: 


and likewise for The ratio between the fc*** and 
[k — I)*** term in the above polynomial is less than 
{af2/a‘ji)i3f{i3f — l)/k = Rf/k. Recalling the definition 
of the Qfn’s, one can roughly estimate it by assuming 
the mass of the hadrons to be the sum of the constituent 
masses of valence quarks and using Eq. 0 : 
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where f = cot b and j / running over all hadrons contain¬ 
ing / or / quark, the partition function can be rewritten 
as: 


Z{Vc, Vc, l^b, Vl) ^ Zl ZeiVe, Vc)Zb{l'b, vi) + E 


3cb 


X Zc{Vc - Vcj, Vc - Vcj)Zb{vb - Vbj,l% - b%j) (5) 


where jcb runs over all hadrons containing both flavours. 
It is now clear from m that we will henceforth be deal¬ 
ing with factorized canonical partition functions because 
both Zc and Zf, do no longer involve hadrons with both 
flavours. Whether stopping the expansion in © at first 
order is sufficient, will be discussed later. If we now de¬ 
note by Ufn the sum of ZjXj for hadrons with n units of 
open flavour /, a/„ for those with n units of open anti¬ 
flavour /, and a/o for // states, we can calculate the Zf 
in © by expanding in powers of a/o: 
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where ges{T,X) is an effective degeneracy parameter in¬ 
cluding the spin degeneracy and the different states with 
the same numbers of heavy quarks, weighted by the ra¬ 
tio ZjXj/ziXi, Zl being the one-particle partition func¬ 
tion of the lowest lying state. In Eq. © we tacitly 
assumed that ^eff is the same for hadrons with one or 
two heavy quarks, which approximately holds according 
to our numerical check. From known states with one 
heavy quark, either c or b, one expects ges 10 at 
T = 165 MeV and Xj = 1. Taking constituent quark 
masses rric = 1.54 GeV, mb = 4.95 GeV, m^^d = 0.33 
GeV ms = 0.51 GeV and T = 165 MeV, i.e. the fit¬ 
ted chemical freeze-out temperature at very large energy 
[1,13, it turns out that Rc ~ — l)/0.34 fm“^V and 

Rb ^ Vbi^b — l)/l-4 fm“^V. Therefore, with the large 
volumes involved in HIC, Rc,Rb are likely to be <C 1 
unless Vc or Ub are consistently large. In this case we 
can approximate the C in Eq. © with its first term, i.e. 
(,( 13 f) ~ a£/i//!, so Eq. (jBJ becomes: 
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Again, the ratio between the and the (/i — 1)*^ term is 
less than a/o^'/^'/Za/iO/j^h, which is approximately equal 
to Rf/h, with Rf quoted in Eq. because 13 f dif¬ 
fers from Vf by few units and being a/i very close to 
afi if chemical potentials are not too large. Therefore, 
under the same conditions needed for its validity, the 
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sum in cnj can be approximated with its first term, i.e. 
Zf{iyf,vjr) ~ ay^ay^/vfliyjrl. By using this approxima¬ 
tion, we can finally estimate the ratio between the first 
and the zeroth order terms in the expansion For 
instance, for mesons carrying one c and one b quark: 
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For hadrons with two b’s and one c and viceversa, it 
can be easily shown that this ratio is even smaller. If 
Rcb 1 (which is likely to be), the first order term in 
the expansion of the partition function 0 is negligible. 
Under these circumstances, and provided that the afore¬ 
mentioned conditions on Rc, Rb are met, the primary 
average multiplicity of heavy flavoured hadrons for fixed 
number of c, c, b, b quarks is especially simple: 

W-A n 


The (I12II is to be further averaged over the multiplicity 
distribution of cc and of bb pairs created in a 
single collision. If they are independently produced, this 
is a Poisson distribution and, for open flavoured hadrons, 
the sum m yields its factorial moments, i.e.: 



= ZjXj 
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whereas for quarkonia it is more complicated since Vc = 
Uc and Ub = vi- Eq. dia is our final formula. As 
has been mentioned, it is an approximated expression 
valid if Rc, Rb, Rcb I. However, it can be shown, 
by performing an asymptotic expansion of the integral 
® [13 that it still holds under the weaker condition 
Rf/vf 1 if Vf is large. Altogether, the Eq. 
implies that the contribution of hadrons carrying more 
than one heavy flavoured quark in the balance equations 
= {vf) is neglected. It is interesting to note 
that the enhancement factors 77 / = {vf)/cLfi are propor¬ 
tional to the density of heavy quarks at the hadroniza- 
tion temperature, as a/i oc V. Therefore, the ratio be¬ 
tween multiply and singly heavy flavoured hadrons, pro¬ 
portional to {{vf), increases with centre-of-mass 
energy because the volume (or the charged multiplicity) 
increases much slower than <Jcc and <jy^ do. 

The formula ca can now be applied to estimate the 
average multiplicity of multiply heavy flavoured hadrons 
in HIC at RHIC and LHC. For sake of simplicity, we 
will confine ourselves to full phase space integrated quan¬ 
tities, disregarding spectra. To get started, we need 
the cross sections tJcc and in pp collisions at rel¬ 
evant energies. There is a large uncertainty on these 


values; recent calculations indicate Ucc = 110 — 656 ph 
and = 1.2 — 2.86 ph at ^/s = 200 GeV and 
= 3.4—9.2 mb and = 88—260 ph at \/s = 5.5 TeV 
[ 1 ^ . The production of heavy quark pairs is a hard pro¬ 
cess and should scale like the number of collisions A^coii 
in the Glauber model. Specifically, if (Jinei is the total 
inelastic NN cross section, the average multiplicity of cc 
pairs is {vf) = {vj) = ff /At RHIC, in Au-Au 
collisions at = 200 GeV, ainei — 42 mb and for a 

5.5% centrality selected sample, corresponding to an im¬ 
pact parameter range 0-3.5 fm, A^coii = 1080 [l^]- Thus, 
the average multiplicity of cc pairs ranges from 2.8 to 
17, whereas for bb pairs from 0.03 to 0.07. On the other 
hand, at LHC, in Pb-Pb collisions at = 5.5 TeV, 

CTinei — 60 mb and for a 5.1% centrality selected sam¬ 
ple, corresponding to the same impact parameter above, 
A^coii = 1670 ^3- l'^ lliis latter case, the average multi¬ 
plicity of cc pairs ranges from 95 to 256 and from 2.4 to 
7.2 for bb pairs. It should be noted that these estimates 
do not take into account possible structure function satu¬ 
ration effects, which are predicted to reduce heavy quark 
cross section at LHC [TJ- Since there is not a clearcut 
evidence of this phenomenon as yet, we stick to the tra¬ 
ditional picture of A^coii scaling, though this possibility is 
worthy of consideration in the future. 

The other key ingredient in our calculation is the vol¬ 
ume V. In order to extrapolate from SPS to RHIC, 
we take advantage of the fact that V is proportional 
to the average multiplicity of pions. In Pb-Pb at SPS 
at = 17.2 GeV, V ~ 3.5 10^ fm^ in full phase 

space |3], and (7r+ -|- 7r“) ~ 1258 ^3- ^1 RHIG, at 
= 200 GeV in full phase space (7r+ -|- 7r“) ~ 3343 
[I 3 I leading to U « 10^ fm^. This extrapolation as¬ 
sumes very little variation of temperature and baryon- 
chemical potentials from SPS to RHIC, which approxi¬ 
mately holds 13 ■ In order to extrapolate from RHIC to 
LHC, we pragmatically use the saturation model (which 
proved to be successful in extrapolating multiplicity from 
SPS to RHIC) which predicts an increase in (uch) by 
a factor ~ 4.5 from = 200 GeV to 5.5 TeV 

[ 13 . This means that Vlhc ~ 4.510"^ fm^ with a 
fair uncertainty up to a factor 2. By using the above 
values for V, and, conservatively, the upper estimates 
for (i/c) and (ub) we obtain, from Eqs. 0 and m 
Rc = 0(10-2), i?f, = O{10-Y,Rcb = 0(10-4) at RHIG 
and Rc = 0(1), Rb = 0(10-4), Rcb = 0(10-4) at LHC, 
by using as input parameters mc,mb,rnu,d and T the 
same quoted below Eq. 0. Therefore, both at RHIC 
and LHC energies the formula 03 should be fairly accu¬ 
rate, because either the R’s are <C 1 or, like in the charm 
sector at LHC, {vc) ^ 1 and Rc/{uc) 1- 

We can now perform our predictions. To estimate 
the 77 /’s, we use the approximation a/i ~ gesV[{mf + 
mu)T/2TrY/^ exp[{—mf + mu)/T] with input parameters 
like for the ratios Rf. We then obtain rjc ~ 1.7—10 
and pb « (3.5 - 8.2) • 10^^ at RHIC and 77 ^ « 12 - 34 
and pb ~ (0.63 — 1.9) • 10® at LHC. With these numbers, 
assuming the mass of multiply heavy flavoured hadrons 
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to be the sum of its quarks constituent masses, using 
Eq. (I Ilf II with Xj = 1 (i.e. taking vanishing chemical po¬ 
tentials) and appropriate spins, we get average primary 
yields of doubly charmed baryons (like Sec and flee) be¬ 
tween 0.7 • 10“"^ and 7 • 10“^ in central collisions at RHIC 
and between 0.019 and 0.38 at LHC. For mixed charmed- 
beautiful hadrons (like Sf,c, f^bc and Be meson), the yields 
should range between 4 • 10“^ and 6 • 10“^ at RHIC and 
3 • 10“^ and 0.022 at LHC. For doubly beautiful baryons 
(like S^b and flbb), our estimates range between 2 • 10“® 
and 3 • 10“® at RHIC and between 2.6 • I0“® and 7 • 10“^ 
at LHC. For the flecc baryon, the predicted yields are 
affected by a large uncertainty due to the cubic depen¬ 
dence on rje] they range between 7 • 10”"^ and 10“"^ at 
RHIC and between 10“^ and 0.03 at LHC. For charmed 
baryon yields at LHC, the predictions of the Eq. 
turn out to be in good agreement with a preliminary 
calculation with the exact formula. All of the previous 
yields are enhanced by the feeding from heavier states, 
even by factor of about 4-5; another factor 2 comes 
from antiparticle yields. These factors should roughly 
compensate for the limited rapidity window accessible 
to experiments. Therefore, while at RHIC only doubly 
charmed hadrons seem to be within reach, at LHC, with 
a statistics of 10^ central events/year, doubly and triply 
charmed, charmed-beautiful, and perhaps doubly beau¬ 
tiful hadrons could in principle be observed. 

We can now compare the above yields with the predic¬ 
tions by production models based on QCD hard scatter¬ 
ing. At LHC, a model where heavy diquarks produced in 
0 {a^) diagrams are assumed to fully hadronize into ccq- 
baryons, yields an upper bound on inclusive production 
of (10-4 - 10-3) (z/c) at Vs = 14 TeV HI in pp, to be 
compared with ges{0.8 — 2) • l{)~^{ve) from coalescence 
in HIC at = 5.5 TeV. A larger difference is found 


in the sector, where in pp at ^/s = 1.8 TeV the IS- 
wave cross-section is predicted to be around 1 nb [l^ . 
implying a ratio {E-be)/(vb) ~ 10-® to be compared with 
(/gg(3—9)-10-4 from coalescence, i.e. at least one order of 
magnitude larger. Since the production process is 0{ag), 
the difference is even larger for fleece for which a recent 
calculation predicts a ratio {flccc)/{vc) = 0(10-^) 
in pp at Vs = 14 TeV; this is between 2-3 orders of 
magnitude lower than our estimated ratio from coales¬ 
cence at Vs^vAT = 5.5 TeV i.e. (0.1 — 1) • IQ-^. For 
charmonia and Be, the two mechanisms give closer pre¬ 
dictions. For {J/tf})/{ve), the difference is estimated to 
be a factor about 2.5 at LHC , whilst for Be cross 
sections calculations |2| at Vs = 14 TeV imply a ratio 
{Be)/{vb) = 0(10-3), around the same as from coales¬ 
cence at VsATAr = 5.5 TeV. Also, it should be pointed 
out that, unlike J/V and Be, multiply heavy flavoured 
hadrons have not been measured in hadronic collisions, 
so the predictions of the models based on hard scattering 
are still to be checked. 

The predominant uncertainty on the previous esti¬ 
mates is that on heavy quark cross section. Other rel¬ 
evant uncertainties are the those on masses, effective de¬ 
generacy, extrapolated temperature and charged multi¬ 
plicities, a modulation of the production as a function of 
rapidity as well as the replacement of the approximated 
formula m with the exact one. Yet, all these effects, 
which will be discussed in detail in a forthcoming paper 
HI , cannot alter the ratios of multiply to singly flavoured 
hadron yields by one order of magnitude. So the conclu¬ 
sion remains that if statistical coalescence scheme applies, 
a large enhancement in the measurement of (Sbc)/(B), 
which becomes dramatic for (flecc) / {D), could be found 
in heavy ion collisions with respect to pp at the LHC en¬ 
ergy. This could be a clear indication of QGP formation. 
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